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Abstract 

Background: Kruppel-like factors (KLFs) are a group of master regulators of gene expression conserved fronn flies to 
human. However, scant information is available on either the mechanisms or functional impact of the coupling of 
KLF proteins to chromatin remodeling machines, a deterministic step in transcriptional regulation. 

Results and discussion: In the current study, we use genome-wide analyses of chromatin immunoprecipitation 
(ChlP-on-Chip) and Affymetrix-based expression profiling to gain insight into how KLFl 1, a human transcription 
factor involved in tumor suppression and metabolic diseases, works by coupling to three co-factor groups: the 
Sin3-histone deacetylase system, WD40-domain containing proteins, and the HPl-histone methyltransferase system. 
Our results reveal that KLFl 1 regulates distinct gene networks involved in metabolism and growth by using single 
or combinatorial coupling events. 

Conclusion: This study, the first of its type for any KLF protein, reveals that interactions with multiple chromatin 
systems are required for the full gene regulatory function of these proteins. 
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Background 

The traditional strategy for studying the role of transcrip- 
tion factors in gene regulation relies on either small or 
large-scale expression analysis following overexpression, 
somatic knockdown, or germ line deletion. While these 
approaches permit the identification of gene networks reg- 
ulated by transcription factors, little information is gath- 
ered about the chromatin coupling events that ultimately 
drive gene expression. This gap in knowledge is unfortu- 
nate as transcription factors are often modular proteins 
armed with multiple sites for potential interactions with 
chromatin cofactors. Thus, knocking-out a transcription 
factor will, in theory, disrupt all chromatin-coupling events. 
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For this reason, it is important to dissect whether different 
transcriptional regulators use single or combinatorial coup- 
ling events to regulate distinct gene expression networks. 

This knowledge is particularly important for proteins 
in which dysfunction contributes to disease. Impairment 
in each of the coupling events may yield differential effects, 
affecting disease penetrance as well as progression. More- 
over, with the development of novel chromatin-centric 
pharmacology, different drugs can partially inactivate 
certain gene networks regulated by a transcription factor 
while leaving others intact. Therefore, the dissection of 
individual and combinatorial chromatin coupling events 
possesses both biological and medical relevance. 

Consequently, the current study focuses on KLF 11 as a 
model system for addressing these important questions. 
KLFll is a well characterized human disease-causing gene 
that couples to several chromatin partners [1]. Alterations 
in KLFll, originally discovered by its role in growth 
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regulation, causes juvenile (MODY7) and neonatal (Ins- 
331 mutation) diabetes [2-4]. KLFll is an inducible gene, 
responsive to a large variety of growth regulatory and 
metabolic stimuli, and functions in the nucleus to regulate 
gene expression by coupling to distinct chromatin part- 
ners. Upon stimuli, KLFll binds to promoters containing 
the consensus CCCCGC/CCCCAC sequences via its three 
C-terminal zinc finger domains and, through well-cha- 
racterized protein-protein interaction modules in its 
N-terminus, differentially recruits chromatin partners 
such as the Sin3-HDAC complex [5], WD40 proteins [6], 
and the HPl-HMT system [7]. Thus, through these 
domains, KLFll translates environmental signals into 
distinct programs of gene expression, which remain to be 
defined in detail. 

In the current study, we have employed a combination 
of genome-wide ChlP-on-Chip and gene expression 
profiles to reconstruct both the direct and indirect ef- 
fects of KLFll on the regulation of different gene 
networks. Furthermore, we used site-directed mutants to 
disrupt individual chromatin coupling domains to gener- 
ate genome-wide expression for the identification of the 
genes regulated by KLFll in response to binding to each 
of its chromatin cofactors. Our results demonstrate that 
certain gene programs require a single chromatin machine 
interaction with KLFll while others require interaction 
with multiple chromatin systems. Accordingly, we provide 
analyses of the distinct gene networks regulated by KLFll 
and its coupling to chromatin partners. This type of 
comprehensive genome-wide analysis has never been per- 
formed for any member of the KLF family. When applied 
to KLFll, our analysis demonstrates that this protein 
behaves in a modular fashion with chromatin machinery 
to alter gene expression and impacts a variety of biological 
processes. 

Results 

Distinct gene expression networks are regulated by the 
differential coupling of KLF1 1 to individual chromatin 
partners 

Recent studies support a model whereby KLFll functions 
by binding to GC-rich sites within promoters of different 
gene networks involved in the regulation of metabolism 
and cell growth [2-4,7-12]. This data is in agreement with 
the biological role of this transcription factor in cancer 
and diabetes. However, how binding of KLFll to these 
promoters regulates gene expression remains to be fully 
understood. Recent data demonstrate that the KLFll pro- 
tein behaves as a scaffold for recruiting different chro- 
matin cofactors via distinct structural motifs (Figure lA). 
The Sin3 Interacting Domain (SID) between amino acids 
22-40 enables coupling of KLFll to the Sin3/HDAC 
system [11-19]. Introduction of proline residues at amino 
acids 29-30 interrupts the association between KLFll 



and the Sin3 scaffold protein, which includes interactions 
with HDACl/2. The region between 281 and 373 amino 
acids is a proline-rich domain that couples with a variety 
of WD40 proteins, including G -protein coupled receptors 
[4]. The A347S mutation observed in MODY7 (neonatal 
diabetes) falls within this domain and has been demon- 
strated to decouple KLFll from novel transcription factor 
Gp2. The zinc finger domain has been demonstrated to 
interact with histone acetyltransferases p300 and CBP 
[3,20,21]. Finally, the region between 483 and 487 dictates 
interaction between KLFll and HP la [7]. Deletion of this 
C-terminal portion decouples HP la from histone methyl- 
transferases SUV39 and G9a, leading to defects in metab- 
olism and tumor suppression. Currently, it is unknown 
whether these systems work in isolation or in a coopera- 
tive manner to regulate gene expression. 

To address this question, we first preformed genome- 
wide promoter binding analysis of wild type KLFll. We 
then limited our analysis to 5 kb upstream and down- 
stream of the transcription start site of known gene pro- 
moters. Examination of the distance of KLF 11 -bound 
peaks from each TSS revealed that the majority (-60%) 
was within 500-1000 base pairs of the TSS (Figure IB). 
Gene ontological analysis of KLFll-bound targets reveals 
association to a large number of metabolic processes, 
including many already characterized as regulated by 
KLFll, such as insulin regulation [3,4,12] and Akt/TOR 
signaling [11] (Figure IC). Next, we performed genome- 
wide expression profiles using mutants for these sites 
which have been shown to specifically uncouple KLFll 
from each of these chromatin remodeling proteins, 
namely EAPP (Sin3/HDAC), A347S (WD40 proteins), and 
A486 (HPla/HMT) (Additional file 1: Figure SI). Of the 
404 gene targets identified by chromatin immunoprecipi- 
tation, for which corresponding gene expression data was 
available, 19% (n = 75) display alteration in the presence of 
the wild type KLFll (p < 0.05). Clustering analysis reveals 
that the expression pattern of these gene targets is 
disrupted in the presence of one or all of the described 
mutants. Thus, we conclude that KLF 11 -mediated gene 
expression may be disrupted by decoupling of the tran- 
scription factor from its chromatin co-factors in a singular 
or combinatorial fashion. 

Distinct gene networks are regulated by the 
combinatorial coupling of KLF1 1 to chromatin partners 

Next, we investigated which gene targets are uniquely 
regulated by each of these pathways. Statistical analyses 
of each expression profile show that, indeed, when com- 
pared with the empty vector control condition, each muta- 
tion induces changes in gene expression that are either 
specific to each chromatin system or similar to the wild 
type. Thus, comparison of the genes significantly affected 
by each condition (p-value <0.05 and fold change ± 1.5) 
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Figure 1 KLF11 -mediated gene expression is disrupted in tlie presence of cliromatin coupling mutations. (A) KLFl 1 is 495 amino acids 
and contains a liiglily conserved C-terminal C2H2 zinc finger domain. Known cliromatin coupling domains explored here are highlighted. 
(B) Panel epithelial cells were transfected with wild type KLFl 1 and chromatin immunoprecipitation performed and hybridized to a whole 
genome promoter array. The frequency of identified peaks 5 kb upstream and downstream of known transcription start sites were binned 
according to distance. The majority of KLFl 1 binding sites were 500-1000 base pairs from transcription start sites. (C) Gene ontological analysis of 
KLFl 1 -bound targets reveals enrichment of genes in known KLFl 1 -associated biological processes, including immune response, TOR signaling, 
and insulin sensitivity. For whole genome analysis, mutants were designed against three of the previously characterized chromatin coupling 
domains. The LAPP mutation in the N-terminus decouples the Sin3/histone deacetylase system. The A347S mutation in the proline rich domain 
decouples KLFl 1 from WD40 containing proteins. Finally, the deletion mutation starting at amino acid 486 disconnects KLFl 1 from the HPl/histone 
methyltransferase system. Panel cells were transduced with empty vector, wild type KLFl 1 or the A347S, A486, or LAPP mutants. Whole genome 
transcriptional profiling was performed. The criteria for significant regulation over empty vector was set at a threshold of +/- 1 .5 log2 fold change and 
a p-value with false discovery rate of less than 0.05. For the LAPP mutant, the p-value did not include false discovery rate thresholding due to the 
limited experimental effects of this mutant. (D) 75 genes were significantly regulated by KLFl 1 overexpression (p < 0.05) and are directly bound by 
KLFl 1 as determined by chromatin immunoprecipitation. Examination of effects of the three chromatin decoupling mutants reveals that expression is 
frequently altered in the presence of one or more of these variants. 



Calvo et a I. BMC Molecular Biology 2014, 15:10 
http://www.bionnedcentral.conn/1471-21 99/1 5/1 0 



Page 4 of 1 5 



shows that 97% of transcripts that are modulated by the 
wild type KLFll are also regulated by any of the three 
mutations, with the A347S and A486 mutations sharing 
the largest number of genes with the wild type (Figure 2A). 
In fact, of the 801 genes modulated by the A486 mutation, 
93.9% are also modulated by the A347S mutation. Finally, 
the three mutants share with KLFll only 11 different 
genes, which likely require an intact coupling of this tran- 
scription factor to all the chromatin proteins involved 
in its function (Figure 2B). We observe that the A347S 
mutation shows the largest number of significantly regu- 
lated transcripts (n = 708) upon the decoupling of KLFll 
from WD40 proteins (Figure 3A). By contrast, the subsets 
of genes uniquely regulated by the A486 or EAPP 
mutants are much smaller, 44 and 21 genes, respectively 
(Figures 3B-C). 

The A486 mutation uniquely modulates only 44 genes 
(5.5%) of the 801 genes significantly modulated under 
expression of the mutant. Almost half of all A486 modu- 
lated genes (n = 390) are common to the A347S mutant 



(Figure 3D) and the other half are common to both 
A347S and wild type KLFll (Figure 3E). Interestingly, 
the A486 and the A347S mutants have, in general, a 
similar direction of modulation, although with varying 
degrees of signal intensity. On the other hand, the 
A347S and A486 modulated genes are almost always 
completely reverted by the EAPP mutation, in which the 
signals are constantly proximal to empty vector values. 

Under the same conditions, the A347S mutant, which 
associates to the development of human juvenile diabetes, 
induces changes in 1521 genes. 46.5% of the modulated 
genes (708/1521) are specific to this mutation, while ap- 
proximately half of the genes (48%) are common to the 
A486 mutant and 43.3% are common to both A347S and 
wild type KLFll (352/423, Figure 3F). Finally, 7% of the 
genes (n = 50) are common to wild type KLFll, and only 
15 common to the EAPP mutant. The EAPP mutation, 
which disrupts coupling to the Sin3-HDAC complex, 
reverses the expression of most genes induced by the wild 
type KLFll or by the other two mutations. 
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Figure 2 KLF11 differentially couples to distinct chromatin 
co-factor systems. (A) Venn diagram of the overlap of genes 
significantly regulated by wild type KLFl 1 and its mutants. Genes 
that occur at the overlap of KLFl 1 and all three of its chromatin- 
coupling mutants were deemed as independent of the effects of 
chromatin de-coupling and therefore inherent to the protein itself. 
(B) Heatmap of the 11 KLFl 1 -regulated genes that are dependent 
on interaction with all three chromatin coupling systems. 



Identification of KLF 11 -regulated genes for which expression 
is independent of known chromatin coupling events 

Subsequently, we examined genes for which expression 
might be regulated by KLFll in a manner that does not 
involve the chromatin remodeling machines examined 
above. A small percentage of genes are still significantly 
modulated uniquely by wild type KLFll (Figure 3G). 
Ontological classification demonstrates that these 14 
genes are associated with diabetes (ALMS IP; RAETIL; 
GYSl; RBM47; EFR3B), cancer (RPL27A; SPAMl; PTCHl; 
ZNF277; USP35; CDC86), or metabolism (C3orfl5; AOC2; 
GOTl). Taken together, these data indicate that the vast 
majority of KLFll s biological processes are dependent on 
its interaction with chromatin regulators within our model 
system. A repertoire of chromatin-independent function is 
predicted, as the transcription factor is capable of inter- 
action with gene promoters directly. Alternatively, the 
regulation of these genes may involve yet unidentified 
chromatin co-factor. 

Pathway reconstruction of combinatorial KLFl 1 -regulated 
chromatin pathways by ontological approach 

Gene Ontology (GO) enrichment analysis of the genes 
modulated by wild type KLFll and each of the three 
mutations show various common biological processes 
(Figure 4A). Indeed, the wild type shares 60% (34/57) of 
the biological processes with at least one of the mutants. 
One of the biological processes, lipid cellular metabolic 
process, is enriched in wild type KLFll and all three 
mutants (Additional file 2: Table SI). Among the most 
highly enriched biological processes, 14 terms are com- 
mon to wild type KLFll and the A347S and A486 muta- 
tions. Variations of the intensity of signal, the number of 



Calvo et ol. BMC Molecular Biology 2014, 15:10 
http://www.bionnedcentral.conn/1471-21 99/1 5/1 0 



Page 5 of 1 5 



A3475 unique genes (n=708) 



B 




A486 unique genes (n=44) ^ EAPP unique genes (n = 21) 





Aa47S - \486 common genes (n-390) 




KLF11-A347S common genes (n=411) 

E 

EAPP 

A48S 
KLF11 
A347S 



KLF11-A347S- \486 common (n=352) 




Figure 3 KLF11 regulates gene expression tiiroughi combinatorial and singular chromatin coupling events. (A) The A347S mutant, which 
decouples KLFl 1 from WD40 proteins, displays the most genome wide effects, uniquely affecting the expression of 708 genes, whereas the (B) 
A486 deletion mutation (decouples KLFll from HPl/histone methyltransferases) and the (C) EAPP mutation (decouples KLFll from binding the 
SinS scaffold protein and subsequently histone deacetylases) uniquely regulate only 44 and 21 genes, respectively. (D) Examination of the overlap 
between genes significantly regulated by the A347S mutant reveals that approximately 50% are regulated in a similar fashion to the A486 
mutant, although with varying degrees of intensity. The EAPP mutant, however, displays near complete reversal of these targets. (E) Wild type 
KLFl 1 and the A347S mutant share 41 1 targets apart from either the EAPP or A486 mutants, although (F) 352 targets exists that are shared 
between all three systems. These data indicate that KLFl 1 chromatin coupling occurs in a largely combinatorial fashion. (G) Only 14 genes were 
identified that are uniquely regulated by KLFl 1, independent of chromatin coupling to the transcription factor. Ontological analysis of these 
genes reveals roles in cancer, cellular proliferation, and metabolism. 



genes implicated, and the score of GO enrichment for 
these 14 biological processes are shown in Additional 
file 2: Table SI. Interestingly, wild type KLFll and 
mutants conserve several biological processes related 
to anabolism or catabolism of lipids, amino acids, and 
glycoconjugates. Other common groups include redox 
processes of protein homodimerization, regulation of 



the activity of protein kinases, and response to the 
estrogen. 

Two of the wild type KLFll biological processes, 
the G2/M transition during DNA damage checkpoint 
and cilium assembly, are also conserved in the A347S 
mutation. There are 16 biological processes exclusively 
conserved between wild type KLFll and the A486 mutant. 
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Figure 4 Biological regulated by KLF11 and its chromatin binding partners. Genes significantly regulated by wild type KLFl 1 and its 
mutants compared to empty vector were analyzed for enrichment of biological processes by an ontological approach. A threshold of 3 genes 
and a p-value of less than 0.05 by Fisher's Exact Test were required to be considered significantly regulated by the transcription factor and its 
mutants. (A) Venn diagram of biological processes demonstrates that the A347S mutant causes the largest number of genome-wide effects with 
significant alterations in 150 processes, 93 of which are unique to the mutant. Only one biological process, lipid cellular metabolism, is dependent 
on the combinatorial effects of all three chromatin coupling systems. The top biological processes were computed for each of the proteins 
compared to empty vector. The pie charts for (B) Wild type KLFl 1, (C) A347S mutant, (D) A486 mutant, and (E) EAPP mutant display the top ten 
scoring biological processes for each with enrichment percentage displayed. These results indicate that the effects of KLFl 1 coupling to chromatin 
co-factors impact on a large number of biological functions. 



These include cellular metabolic processes, such as alde- 
hyde and nitrogen compound metabolism, isocitrate and 
2-oxoglutarate metabolism, glycosaminoglycan catabolic 
processes, and fatty acid beta-oxidation. 



Interestingly, the A347S and A486 mutations share 
more than a quarter of biological processes (n = 39) that 
are absent in wild type KLFll and the EAPP mutant. 
Between these 39 GO terms, we found enrichment of 
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terms associated to the monitoring of the DNA transcrip- 
tion, RNA polymerase I transcription promoter regulation, 
and epigenetic control of gene expression, such as histone 
H3 acetylation, chromatin remodeling, protein acetylation 
and deacetylation, and regulation of phosphorylation. 
Wild type KLFll and EAPP share only two GO-terms in 
common: branched chain amino acid family catabolic 
processes and negative regulation of TOR and signaling 
cascade. 

Overall, 22 of 57 (38.6%) of the biological processes 
enriched in the wild type KLFll are not present any of 
the three mutants. Among the 22 most enriched bio- 
logical process found only in the wild type KLFll, we 
observed several metabolic processes, such as glucose, 
carbohydrate, and very long-chain fatty acid metabolism, 
oxaloacetate metabolism, fatty acid homeostasis, and 
oligosaccharide metabolism. Top processes associated 
with wild type KLFll expression are shown in Figure 4B 
and Additional file 3: Table S2. Other specific mechanisms 
of wild type KLFll lost by mutation are two mechanisms 
associated with gene expression (transcription from the 
RNA polymerase III promoter and mRNA capping) and 
two biological processes associated with the integrity of 
the DNA (the DNA repair and the response to DNA dam- 
age stimulus). Finally, three important biological processes 
(e.g., the mitotic cell cycle, the negative regulation of 
insulin-like growth factor receptor signaling pathway, and 
the negative regulation of epithelial cell proliferation) are 
also lost by the mutants. From these data, we conclude 
that the combinatorial association between the transcrip- 
tion factor and a variety of chromatin coupling systems 
dictates the majority of KLFll function, particularly in 
regulating gene expression and metabolism. The high de- 
gree of overlap suggests a form of "epigenetic redundancy" 
to ensure the proper regulation of these gene targets and 
processes. 

Pathway reconstruction of singular KLF1 1 -regulated 
chromatin pathways by ontological approach 

Of the mutants, A347S displays the highest number of 
unique biological process that are not enriched in the 
wild type (n = 93). Top processes associated with A347S 
mutant expression are shown in Figure 4C and Additional 
file 3: Table S2. Interestingly, almost 30% of these biological 
processes are associated with the structure and mainten- 
ance of DNA or telomeres (Additional file 4: Table S3, 
List 1), processes related to phosphorylation and dephos- 
phorylation, deacetylation, methylation, and demethylation 
(Additional file 4: Table S3, List 2), and processes of 
transcriptional regulation by poly-II and mRNA transport 
(Additional file 4: Table S3, List 3). 

In turn, the A486 mutant exclusively enriches 45 
biological processes that are mediated by the wild type 
KLFll or in other mutants. Top processes associated 



with A486 expression are shown in Figure 4C Additional 
file 3: Table S2. Unlike the A347S mutant, mutation 
A486 is particularly enriched in processes associated 
with metabolic and biosynthetic processes, protein ex- 
port-import, and energy protein modification associated 
processes (Additional file 5: Table S4, List 1). Other 
groups of processes specifically enriched in this mutant 
are associated with proliferation and epithelial cell survival 
(Additional file 5: Table S4, List 2). 

The EAPP mutant did not display robust gene expres- 
sion changes within the experiment. However, utilizing a 
more permissive cutoff yields a list of 53 biological cat- 
egories, although each only contains a single gene. We 
find enrichment of biological processes involving nitric 
oxide transport, catabolic processes, and developmental 
morphogenesis function, among others. Top processes 
associated with EAPP expression are shown in Figure 4E 
and Additional file 3: Table S2. From these data, we are 
able to identify singular KLFll chromatin coupling events 
dictating a number of critical processes, particularly in 
management of cell cycle control and DNA replication, 
ones that inherently require a more specific and finely 
regulated degree of control to execute their programs 
under strict spatial and temporal constraints. 

Identification of canonical signaling pathways mediated 
by KLF1 1 and mutants 

Subsequently, we wished to further assess the relatedness 
between the gene targets significantly modulated by KLFll 
and its mutants utilizing a semantic-based approach and 
the Ingenuity Global Canonical Pathways (GPC) algorithm. 
The input dataset of significantly regulated gene targets 
was compared against the IPA canonical pathways, which 
are curated from published literature. The significance of 
the associations between the data set and a given canonical 
pathway is determined by the ratio of the number of genes 
mapping to the pathway divided by the total number of 
pathway genes. A p-value is calculated using Fischer's 
Exact Test determining the probability that the association 
between the data set and the pathway occurs by chance 
alone. GPC analysis of genes modulated by the wild type 
KLFll and the three mutants shows a variety of signifi- 
cantly enriched pathways (p-value <0.05). The wild type 
KLFll shares about 71% (30/42) of its pathways with at 
least one of the mutants (Figure 5A). The most important 
intersection of common pathways is observed between 
wild type KLFll and mutants A347S and A486. Indeed, 
90% (38/42) of these pathways are also enriched in at least 
one of these two mutants. Only two pathways are shared 
between the wild type and the EAPP mutant. 

Nevertheless, all genes implicated in these pathways 
are differently modulated by wild type KLFll and each 
one of the mutants. The expression of these genes signifi- 
cantly regulated by wild type KLFll and the A347S and 
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Figure 5 Canonical signaling pathways mediated by KLF11 and mutants. Using the Ingenuity Global Canonical Pathways algorithm, 
significantly modulated signaling pathways were identified for wild type KLFl 1 and its mutants. A p-value of less than 0.05 as determined by 
Fisher's Exact Test was used as criteria for significant association of focus genes to predetermined pathways curated from published literature. 
(A) Venn diagram of significant canonical pathways regulated by wild type KLFl 1 and mutants reveals that the effects of the EAPP mutation are 
less significant in this experiment compared to the A347S and 486 mutants. These data are supported by the clustering of all significantly 
regulated targets as shown in (B) which demonstrated that the A347S and A486 mutants cluster closely with wild type KLFl 1 while the EAPP 
mutant clusters with the empty vector control, indicating complete reversal of the effects of the other proteins. (C) Among the 28 top scoring 
unique canonical pathways mediated by the 374 mutant are signaling pathways, including, EGF, PDGR, p53, RANK and reelin signaling, and 
glucocorticoid receptor signaling. The A486 (D) and EAPP (E) mutants mediate 4 and 6 pathways respectively, both of which are centered on 
processes involved on the degradation of biogenic amines. These data support KLFl 1 as a master transcriptional regulator that mediates large, 
interconnected signaling cascades. 



A486 mutants are almost always modulated in the 
same sense and clustered together (Figure 5B). In 
contrast, the EAPP site produces a complete reversion 
of the signals with the levels of gene expression close 
to the empty vector level, as evidence by its clustering 



with empty vector apart from the wild type KLFll and 
other mutants. 

Twelve of the 42 significantly enriched pathways are 
conserved only in the wild type KLFll. These normal 
functions of the wild type KLFll are potentially lost by 
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the three studied mutations. The most enriched of these 
pathways includes the nucleotide excision repair path- 
way and the alpha-adrenergic signaling pathways. Inter- 
estingly, of all the remaining pathways, 80% (9/11) are 
metabolic biosynthetic pathways. 

Mutant A347S possesses the most exclusive pathways. 
Indeed, 28 pathways are enriched in the A347S mu- 
tant compared to wild type KLFll and other mutants 
(Figure 5C). Interestingly, over 30% (9/28) of these path- 
ways are signaling pathways, including, EGF, PDGR, p53, 
RANK and reelin signaling, and glucocorticoid receptor 
signaling. Furthermore, pathways involved in diabetes and 
pancreatic adenocarcinoma are also present. These path- 
ways may explain the massive number of genes differen- 
tially modulated by this mutation (n = 1521 genes) and 
suggests that this mutant could induce new potentially 
unexpected and significant physiopathological changes. 

86% of the enriched pathways in the A486 mutant are 
shared with the wild type KLFll or the other mutants. 
Only four new pathways are specifically enriched by the 
A486 mutation (Figure 5D). Interestingly, these pathways 
are implicated in the degradation of biogenic amines 
(an important field associated to physiopathology of 
depression and metabolism of drugs of abuse). Finally, 
the EAPP mutant influences only six pathways exclusively 
(Figure 5E), also primarily related to the degradation 
of biogenic amines. Taken together, we conclude that 
KLFll chromatin-coupled gene targets delineate into 
specific, well-ordered pathways, consistent with KLFlls 
position as a master transcriptional regulator that me- 
diates large cascades of other transcription factors 
and regulators. The effects of KLFll chromatin coup- 
ling are therefore not incidental to its function but a 
primary means by which to execute gene activation or 
repression. 

Reconstruction of downstream biological and disease 
networks mediated by KLF1 1 and mutants 

We next explored the degree of connectivity between 
KLFll and mutant-regulated genes and their pathobio- 
logical associations. Networks of significantly regulated 
genes were algorithmically generated based on their 
connectivity and assigned a score that encapsulates the 
relevance of the generated network to the original list of 
focus genes. All edges are supported by at least one 
literature reference of direct physical, transcriptional, 
and enzymatic interactions, or from canonical informa- 
tion stored in the Ingenuity Pathways Knowledge Base. 
A right-tailed Fisher s Exact Test was used to calculate 
the p-value for networks and a threshold of p < 0.05 used 
to determine significance. A functional analysis of a net- 
work then determined the biological functions and/or 
diseases that are most significant to the genes contained 
within the network (Figure 6). 



Of the identified processes, only three were common 
to wild type KLFll and all three mutants: cell death and 
survival, hereditary disorders, and molecular transport 
mechanisms. Wild type KLFll possessed the most unique 
networks (n = 10) and included a number of characterized 
biological and disease networks, including cellular growth 
and proliferation [22], endocrine disorders, reproductive 
system development and function [23], and biliary hy- 
perplasia, among others. The A347S mutant uniquely 
associates to 6 networks, including ones linked to auditory 
disease, cellular compromise, and organismal development. 
The A486 mutant uniquely associated to 8 networks. 
Among the A486 mediated networks are significant 
associations to gastrointestinal disease, connective tissue 
development, and immune presentation and response. 
Together, these data recapitulate the repertoire of known 
KLFll-mediated functions and uncover a number of new 
potential functional and disease associations for future 
studies. 

Identification of upstream regulators of KLF11 -chromatin 
coupled pathways 

Finally, we wished to identify the upstream transcrip- 
tional regulators using Ingenuity's Upstream Regulator 
(UR) analysis. The UR analytic is based on prior know- 
ledge of the expected effects of transcriptional regulators 
and their target genes stored in the Ingenuity Public 
Knowledge Base. The algorithm examines how many 
known targets of each transcriptional regulator are present 
in the list of significantly regulated targets and compares 
direction of expression to that expected from previously 
published data. For each potential regulator, an overlap 
p-value using Fishers Exact Test and an activation z-score 
was computed, the latter of which determines the activa- 
tion state. A threshold p-value of <0.01 and activation 
score of +/- 2 was considered significant for our purposes. 
The end result was the identification of upstream regula- 
tors of KLFll or its mutants that permit the generation of 
plausible signaling cascades mediated by these upstream 
regulators through KLFll chromatin coupling. 

The results of the UR analysis are presented in 
Additional file 6: Table S5. Wild type KLFll displayed 
two significant upstream inhibitory regulators, PPARG 
and HNFIA, both of which have previously been impli- 
cated in the regulation of the transcription factor or as a 
co-regulator (Figure 7 A) [24]. The A347S mutant, con- 
sistent with the gene level, pathway, and network analyses, 
is the most disruptive, possessing 9 activating upstream 
regulators and 23 inhibitory upstream regulators. Interest- 
ingly, EZH2, a histone methyltransferases, is identified as 
potential upstream activator. Previous reports have identi- 
fied an antagonistic relationship between KLFIO and the 
functions of EZH2/Polycomb Repressive Complex [21]. 
This data hints to other important chromatin couplings to 
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Figure 6 Biological disease networks mediated by KLF11 and its mutants. Using Ingenuity's networl< analysis tool, significantly regulated 
genes by KLFl 1 (A) and its A347S (B) and A486 (C) mutants were examined for connectivity. Networks were scored by relevance and a p-value 
assigned by Fisher's Exact Test to discriminate networks generated by chance alone. A threshold of p < 0.05 was employed to determine signifi- 
cant networks. Networks were then assessed for association to known biological or disease processes. Wild type KLFl 1 and its mutants share three 
processes: death and survival, hereditary disorders, and molecular transport. Of the identified processes, 10 are unique to wild type KLFl 1 and 6 and 8 are 
unique to the A347S and A486 mutants, respectively. A number of known KLFl 1 -mediated biological and diseases processes were revealed by this analysis, 
providing an internal control, including endocrine disorders, gastrointestinal disorders, and cellular growth and proliferation processes. However, the analysis 
also generated a number of novel biological and disease processes that remain to be experimentally validated for KLFl 1. 



KLFll that remain to be completely elucidated. Figure 7B 
shows the genes affected in the presence of the mutant 
when 486 is inhibited by SREBFl or JUN and ATF4. For 
the 486 mutant, 6 upstream activators were identified and 
15 upstream inhibitors, including SREBF2 and HIFIA 
(Figure 7C). Although the EAPP did not display any 



significant upstream regulators as defined by our full 
criteria, a number of other KLF proteins (KLFl, 3, 7, 13) 
were identified as upstream regulators by p-value alone, 
suggesting that inter-association of KLF proteins is also 
critical to KLFll function. In addition, we find that KLFIO 
is identified as affected downstream of both the 486 and 
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Figure 7 Upstream regulators of KLF11 and its mutants. Gene targets significantly regulated by KLFl 1 and its mutants were analyzed by 
Ingenuity's Upstream Regulator Analytic that compares the experimentally derive activation or inhibition of focus molecules to relationships 
between upstream regulators and target molecules known from published data. For each potential upstream regulator and its targets, a p-value 
for the degree of overlap between that of the protein under study and a z-score for the activation status are calculated. For our study, a threshold 
of p < 0.05 and a minimum z-score of +/- 2 were employed. The full results are published in Additional file 6: Table S5 and sample signaling 
cascades presented for wild type KLFl 1 (A), A347S (B), and A486 (C). No cascades were generated that proved significant for the EAPP mutant 
although a more permissive criteria of a p-value of less than 0.05 did generate a number of single gene associations listed in Additional file 6: 
Table S5, including a variety of other KLF proteins, hinting at co-regulation within the KLF family. Many of the upstream regulators, such as PPARG, 
are previously identified co-regulators of KLFl 1. While the upstream regulators frequently overlap between wild type and mutant conditions, the 
end gene targets are unique to each condition. These results suggest that the deregulation of KLFl 1 from a large number of gene networks is 
completely dependent on the effective coupling of the transcription factor to its chromatin co-factor system. 



A347S mutants through multiple KLFll upstream regula- 
tors. Together, these data suggest that transregulation 
among KLF proteins is critical to their function and reveal 
the potential for a chromatin-dependent interplay between 
family members, permitting additional specificity in the 
regulation of distinct nodes within their targeted gene ex- 
pression networks. Although a number of the regulators 
are common to the wild type KLFll and the A347S and 
A486 mutants, indicating that mutation does not neces- 
sarily impair upstream regulation or co-regulation, the ac- 
tivation or inhibition of these proteins regulates a distinct 
subset of endpoint genes. The absence of wild type KLFll 
targets and acquisition of new, inappropriate targets in the 
presence of mutation places KLFll chromatin coupling 
events at the center of the transcription factor s function 
or dysfunction. 

In summary, here we have identified the precise inter- 
play between KLFll and three of its chromatin coupling 
systems, namely Sin3/HDAC, HPl/HMTs, and WD40 
proteins, on a genome-wide scale. The functional repertoire 
of KLFll appears to be dictated by both combinatorial 
effects between these systems and singular chromatin coup- 
ling events. Redundancy, particularly in critical biological 



processes such as metabolism and organismal develop- 
ment, occurs between systems to ensure preservation of 
critical functions. On the other hand, singular chromatin 
coupling events exist to ensure precision in processes 
requiring strict spatial and temporal parameters to avoid 
inappropriate signaling. 

Discussion 

Until recently, the relationship between chromatin and 
transcription factors appeared largely passive. However, 
the evidence that overexpression of four transcription 
factors— Oct4, Sox2, m-Myc, and KLF4— is sufficient to 
reprogram fully differentiated fibroblasts into pluripotent 
ES-like cells suggests an active role for transcription 
factors in the epigenetic regulation of chromatin state 
[25]. KLF proteins, in particular, have demonstrated a 
multitude of feed-forward effects on epigenetic regulation 
of gene expression. Here, we propose KLFll as a model 
for understanding the critical role of sequence-specific 
factors in epigenetic regulation of a multitude of pathobio- 
logical processes. 

KLF proteins deliver epigenetic information to gene 
promoters through three primary mechanisms: (i) sensing 
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and translating environmental stimuli into a program of 
gene expression, (ii) sequence-specific targeting of chro- 
matin remodeling complexes to gene promoters, and (iii) 
transactivation of other transcription factors to assist in 
the regulation of large networks of interdependent genes. 
Transient regulation of gene expression occurs through 
association with histone acetyltransferases (HATs) and 
histone deacetyltransferases (HDACs), while long-term 
gene silencing is enacted by interactions with histone 
methyltransferases (HMTs) and DNA methyltransferases 
(DNMTs). No single classification scheme may accurately 
characterize the function of KLFll in its entirety. How- 
ever, the classification of family members by virtue of its 
co-factors provides a framework to understand the func- 
tional differences in the manner in which it engages chro- 
matin to activate or repress transcription in a dynamic 
and reversible manner. 

The association between KLFll and chromatin coupling 
is well established. Initial biochemical characterization of 
KLFll revealed that its N- terminus possesses a domain 
that adopts an a-helical conformation and that mutations 
within this domain significantly disrupted its binding to 
Sin3 and the subsequent recruitment of HDACs [26,27]. 
Subsequently, KLFll was found to associate with HPla 
[7], one of the epigenetic "gatekeepers" of gene silencing 
through its extreme C-terminus [28-30]. HPl proteins 
repress gene expression by binding to H3K9me marks and 
interacting to H3K9 histone methyltransferases, such as 
G9a or SUV39H1, which methylate this same residue on 
adjacent nucleosomes [31,32]. Deletion of the C-terminus 
leads to deregulation of tumor suppression functions 
mediated by KLFll. Finally, KLFll has also been demon- 
strated to associate with WD40, WWI, WWII, and 
SH3-domain containing proteins through a proline- 
rich domain [4]. Interaction of KLFll and WD40 protein 
G|32, for example, is disrupted in the presence of the 
A347S mutation, the variant associated with development 
of Maturity Onset Diabetes of the Young 7 (MODY7). 
The effect of decoupling of KLFll to chromatin co- 
factors, however, has only been studied on small subsets 
of genes related to the system of focus. 

The investigation presented here is the first to attempt 
to dissect the roles of KLFll and chromatin coupling on 
a genome-wide level to ascertain the combinatorial or 
unique effects of each type of previously identified inter- 
acting chromatin system. Utilizing a single cell type with 
over expression of wild type KLFll and the three previ- 
ously described mutants, A347S, A486, and EAPP, we 
were able to examine the interplay between coupling of 
the transcription factor to the WD40, HPl/HMT, and 
Sin3/HDAC systems, respectively. By introducing saturat- 
ing amounts of each mutant into pancreatic epithelial cells, 
we were able to tease apart the relative and overlapping 
contributions of each chromatin coupling mechanism to 



the regulation of pathobiological gene networks. Our 
results reveal that in the pancreatic epithelial cell utilized 
as our model, a large number of the affected genes, 
networks, and signaling cascades are modulated by the 
coupling of KLFll to its variety of chromatin systems. The 
number of genes identified as regulated independently of 
these systems represented <1% of the total genes mediated 
in the presence of wild type or mutant proteins. 

It is important to underscore that the study reported 
here is unique as it uses genome-wide approaches to dis- 
sect the function of several different chromatin-based 
pathways in the transcriptional function of a single KLF 
protein. The KLF family of transcription factors consists 
of 17 members [1]. Sequence identity at the carboxyl 
terminus among KLF family members is greater than >65%, 
suggesting the regulation of similar types of gene promoters 
and the potential for synergistic or antagonistic regulation 
by family members. Despite the high degree of sequence 
identity in their DNA-binding ability, the functional activ- 
ities of KLF members differ widely due to the high degree 
of variability within the N-terminal domain of each pro- 
tein, the region responsible for interaction with chromatin 
machinery. Although the highly conserved DNA-binding 
domain denotes their familial origins, the ultimate func- 
tional identity of each KLF protein rests in its their ability 
to couple to the molecular machinery that regulate gene 
expression, in particular to bind writers, readers, and 
erasers of the histone code. Although some family mem- 
bers bind to the same chromatin partners as KLFll [33], 
the degree to which the gene networks and biological 
processes identified here may be generalized to other 
family members remains to be elucidated. However, it is 
tempting to speculate that the most closely related pro- 
teins KLF9, KLFIO, KLFll, KLF13, KLF14, and KLF16 
being extensively studied in our laboratory because of 
their participation in metabolic processes, may at least 
regulate some overlapping gene networks. These proteins 
contain distinct motifs which couple to similar and also 
distinct chromatin remodeling complexes. In fact, it has 
been shown that these proteins couple to the same 
chromatin regulators in organisms ranging from flies 
to human. Notably, disruption of the KLFll orthologous 
fly gene Cabut, leads to metabolic problems [34]. In 
human, alterations in KLFll causes neonatal and juvenile 
Diabetes while its closely related gene, KLF14 causes 
obesity and adult type II diabetes [35]. Interestingly, this 
function of KLF 14 appears to be influenced by a cross 
regulation with KLF13. This is interesting since our 
microarray data shows that, KLF9 is significantly down- 
regulated by wild type KLFll and KLFIO is upregulated 
by the 486 and A347S mutants. Therefore, this data sug- 
gests that this subgroup of highly related subfamily of pro- 
teins, which appears to have evolved from the fly ancestor 
Cabut, share similar chromatin coupling mechanisms. 
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similar general function in metabolism, and may, to some 
extent, share targeted gene networks However, they still 
have some degree of biochemical differences as well as 
phenotypic outcomes. On the other hand, the rest of the 
KLF proteins, which are structurally and functionally less 
related to KLF 11 but may still share some chromatin part- 
ners, are predicted to have fewer overlapping co-regulated 
gene networks. In this regard, we are optimistic that the 
type of investigations reported here will fuel the interest of 
both basic biologists and translational scientists who are 
keen on better understanding the function of this exten- 
ded family of transcription factor proteins as it relates to 
biochemistry, genetics, epigenetics, general biology and 
the pathobiology of human diseases. 

Conclusion 

More than fifteen years ago, our research team hypothe- 
sized that the discovery and study of KLF proteins and 
their chromatin cofactors would assist in unraveling 
complex human diseases. In the intervening years, our 
understanding of the interplay between KLF proteins and 
epigenetic machinery in transient and long-term gene 
regulation has grown exponentially. With knowledge in- 
ferred from the predictive power of rationally derived 
computational models, we propose a new paradigm for 
KLF regulation of gene networks through the translation 
of input from cellular milieu into epigenetic information 
to affect changes in chromatin structure. KLF proteins 
currently provide the best model for understanding the 
interactions between mechanisms of "hard inheritance" 
(environmental and genetic variation) and "soft inherit- 
ance" (epigenetic variation) in underscoring the phenotypic 
variability observed in complex disease mechanisms. Per- 
haps most importantly, these new pathways open up infin- 
ite possibilities for targeting new molecules and processes 
in therapeutic intervention and disease management. 

Methods 

Cell culture 

Panel cell lines were obtained from the American Type 
Culture Collection (ATCC, Rockville, MD). Cells were 
cultured as described previously [4,36]. 

Constructs 

Standard molecular biology techniques were used to 
clone full length KLFll, KLF11-A347S, KLFll-EAPP, 
and KLFll-486 into pcDNA3.1/His (Invitrogen, Carlsbad, 
CA). All constructs were verified by sequencing at the 
Mayo Clinic Molecular Biology Core Facility. Epitope- 
tagged (6XHis-Xpress™) KLFll, KLF11-A347S, KLFll- 
EAPP, KLFll A486 variants as well as empty vector 
(AdSCMV) were generated as recombinant adenovirus in 
collaboration with the Gene Transfer Vector Core at the 
University of Iowa. 



Genome-wide expression profiles of KLF1 1 and mutants 

Panel epithelial cells were plated at a density of 1x10^ 
cells/ 100 mm dish and transduced with empty vector, 
KLFll, KLF11-A347S, KLFll- A486, or KLFll-EAPP at 
an MOI of 150. Relative expression levels were con- 
firmed by Western Blot (Additional file 1: Figure SI). 
RNA was prepared as previously described from pooled 
biological triplicates [37]. Global gene expression profil- 
ing was carried out in technical duplicate at the Micro- 
arrays Facility of the Research Center of Laval University 
CRCHUL utilizing the Aff)^metrix Human Gene 1.0 ST 
arrays (28,869 well-annotated genes and 764,885 distinct 
probes). Intensity files were generated by Affymetrix 
GCS 3000 7G and the Gene-Chip Operating Software 
(Aff)^metrix, Santa Clara, CA). A subset of genes was 
validated by qPCR as previously described (Additional 
file 7: Figure S2) [20,38]. 

Genome-wide promoter binding profile of KLF1 1 

Panel epithelial cells were transfected with full-length 
His-tagged KLFll. ChIP was performed as previously 
described [38-41] using an antibody against the His-Tag 
(OMNI D8; Santa Cruz Biotechnology) to detect recom- 
binant expression of KLFll. Non-specific IgG antibody 
was utilized as a negative control. Binding activity was 
derived using the NimbleGen human promoter hybridi- 
zation system (Madison, WI). Peaks were detected by 
searching for >4 probes where signals were above the 
specified cutoff values (90% to 15%) using a 500 bp sliding 
window along 5 kb upstream and downstream of the tran- 
scriptional start site in human promoters. Each peak was 
assigned a score that is the log2 ratio of the fourth highest 
probe in each peak. If multiple peaks are present, the peak 
nearest the TSS is reported. Ratio data was then random- 
ized 20 times to evaluate the false discovery rate (FDR). 
Only peaks with FDR scores <0.2 were deemed high confi- 
dence binding sites and reported. 

Data analysis 

Data analysis, background subtraction and intensity 
normalization was performed using Robust Multiarray 
Analysis (RMA) [42]. Genes that were differentially ex- 
pressed along with false discovery rate (FDR) were 
estimated from t test (>0.005) and corrected using Bayes 
approach [43,44]. A threshold of +/- 1.5 log2 fold change 
with a p-value with FDR of less than 0.05 (without FDR 
for LAPP mutant) was used to determine significantly 
regulated targets. Data analysis, hierarchical clustering, 
and ontology were performed with the OneChanelGUI to 
extend affylmGUI graphical interface capabilities and 
Partek Genomics Suite, version 6.6 (Partek Inc., St. Louis, 
MO) with ANOVA and GO ontological analysis [45]. 
Selected probes and their fold changes were loaded into 
Ingenuity Pathways Analysis Software (Ingenuity® Systems, 
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www.ingenuity.com) for annotation, redundancy checks, 
canonical pathway, biological network, and upstream 
regulator analysis using default parameters. 

Availability of supporting data 

The data sets supporting the results of this article are 
available in the NCBI repository, GSE56778. 

Additional files 



Additional file 1: Figure SI. Analysis of wild type KLFl 1 and mutant 
protein expression. Protein expression of epitope-tagged (Hiss-Xpress™) 
wild type KLFl 1 and the A347S, A486, and EAPP mutants demonstrating 
similar overexpression in Panel cells. a-Tubulin was used as a loading 
control. 

Additional file 2: Table SI. Overlapping biological function mediated 
by KLFl 1 and mutants. 

Additional file 3: Table S2. Top biological processes mediated by 
KLFl 1 and mutants compared to empty vector. 

Additional file 4: Table S3. Biological processes unique to the 
decoupling of KLFl 1 from WD40 proteins (A347S mutant). 

Additional file 5: Table S4. Biological processes unique to the 
decoupling of KLFl 1 from HPla/HMT (A486 mutant). 

Additional file 6: Table S5. Upstream regulators of KLFl 1 and mutants. 

Additional file 7: Figure S2. qPCR validation of Affymetrix gene 
expression data. A small subset of significantly regulated genes identified 
by Affymetrix whole-genome microarray for wild type KLFl 1 or the 
A347S, A486, and EAPP mutants were validated by qPCR. 



Competing interests 

The authors declare that no competing interests exist 
Authors' contributions 

EC, AG, GL, and AM carried out the molecular, promoter binding, and gene 
expression studies, participated in the bioinformatics analysis, and drafted 
the manuscript. Jl and RU conceived of the study, participated in its design 
and coordination, and helped to draft the manuscript. All authors read and 
approved the final manuscript. 

Acknowledgements 

This work was supported, in whole or in part, by National Institutes of Health 
Grants DK52913 (to R. U.), T32-CA1 48073 (to A. G.), Mayo Clinic Center for 
Cell Signaling in Gastroenterology (P30DK084567) and a Career Development 
Award from the Mayo Clinic SPORE in Pancreatic Cancer (P50 CAl 02701) (to G.L). 

Author details 

^Molecular Endocrinology and Oncology Research Center, CHUL Research 
Center, Quebec, Canada. ^Laboratory of Epigenetics and Chromatin 
Dynamics, Mayo Clinic, Rochester, MN 55905, USA. ^INSERM U.624, Stress 
Cellulaire, 163 Avenue de Luminy, Case 915, Pare Scientifique et 
Technologique de Luminy, 13288 Marseille, Cedex 9, France. ^Translational 
Epigenomics Program, Center for Individualized Medicine (CIM), Mayo Clinic, 
Rochester, MN 55905, USA. ^Departments of Medicine, Physiology and 
Biochemistry, Mayo Clinic, 200 First Street SW, Guggenheim 10, Rochester, 
MN 55905, USA. 

Received: 23 September 2013 Accepted: 7 May 2014 
Published: 25 May 2014 

References 

1. Lomberk G, Urrutia R: The family feud: turning off Spl by Spl-like KLF 

proteins. Biochem J 2005, 392(Pt 1):1-1 1. 

2. Bonnefond A, Lomberk G, Buttar N, Busiah K, Vaillant E, Lobbens S, Yengo L, 
Dechaume A, Mignot B, Simon A, Scharfmann R, Neve B, Tanyolac S, 
Hodoglugil U, Pattou F, Cave H, lovanna J, Stein R, Polak M, Vaxillaire M, 



Froguel P, Urrutia R: Disruption of a novel Kruppel-like transcription factor 
p300-regulated pathway for insulin biosynthesis revealed by studies of 
the C.-331 INS mutation found in neonatal diabetes mellitus. J Biol Chem 
2011, 286(3 2):2841 4-28424. 

3. Fernandez-Zapico ME, van Velkinburgh JC, Gutierrez-Aguilar R, Neve B, 
Froguel P, Urrutia R, Stein R: M0DY7 gene, KLF11, is a novel 
p300-dependent regulator of Pdx-1 (M0DY4) transcription in 
pancreatic islet beta cells. J Biol Chem 2009, 284(52):36482-36490. 

4. Lomberk G, Grzenda A, Mathison A, Escande C, Zhang JS, Calvo E, Miller LJ, 
lovanna J, Chini EN, Fernandez-Zapico ME, Urrutia R: Kruppel-like factor 11 
regulates the expression of metabolic genes via an evolutionarily 
conserved protein-interaction domain functionally disrupted in maturity 
onset diabetes of the young. J Biol Chem 2013, 288(24):1 7745-1 7758. 

5. Grzenda A, Lomberk G, Zhang JS, Urrutia R: Sin3: master scaffold and 
transcriptional corepressor. Biochim Biophys Acta 2009, 1 789(6-8):443-450. 

6. Grzenda A, Lomberk G, Svingen P, Mathison A, Calvo E, lovanna J, Xiong Y, 
Faubion W, Urrutia R: Functional characterization of EZH2beta reveals the 
increased complexity of EZH2 isoforms involved in the regulation of 
mammalian gene expression. Epigenetics Chromatin 2013, 6(1 ):3. 

7. Lomberk G, Mathison AJ, Grzenda A, Seo S, DeMars CJ, Rizvi S, 
Bonilla-Velez J, Calvo E, Fernandez-Zapico ME, lovanna J, Buttar NS, Urrutia 
R: Sequence-specific recruitment of heterochromatin protein 1 via 
interaction with Kruppel-like factor 11, a human transcription factor 
involved in tumor suppression and metabolic diseases. J Biol Chem 2012, 
287(1 6):1 3026-1 3039. 

8. Seo S, Lomberk G, Mathison A, Buttar N, Podratz J, Calvo E, lovanna J, 
Brimijoin S, Windebank A, Urrutia R: Kruppel-like factor 11 differentially 
couples to histone acetyltransferase and histone methyltransferase 
chromatin remodeling pathways to transcriptionally regulate dopamine 
D2 receptor in neuronal cells. J Biol Chem 2012, 287(1 6):1 2723-1 2735. 

9. Fan Y, Guo Y, Zhang J, Subramaniam M, Song CZ, Urrutia R, Chen YE: 
Kruppel-like factor-1 1, a transcription factor involved in diabetes 
mellitus, suppresses endothelial cell activation via the nuclear 
factor-kappaB signaling pathway. Arterioscler Thromb Vase Biol 2012, 
32(1 2):2981 -2988. 

10. Fernandez-Zapico ME, Lomberk GA, Tsuji S, DeMars CJ, Bardsley MR, Lin YH, 
Almada LL, Han JJ, Mukhopadhyay D, Ordog T Buttar NS, Urrutia R: A 
functional family-wide screening of SP/KLF proteins identifies a subset 
of suppressors of KRAS-mediated cell growth. Biochem J 20]], 
435(2):529-537. 

1 1 . Buttar NS, DeMars CJ, Lomberk G, Rizvi S, Bonilla-Velez J, Achra S, Rashtak S, 
Wang KK, Fernandez-Zapico ME, Urrutia R: Distinct role of Kruppel-like 
factor 1 1 in the regulation of prostaglandin E2 biosynthesis. J Biol Chem 
2010, 285(1 5):1 1433-1 1444. 

12. Neve B, Fernandez-Zapico ME, Ashkenazi-Katalan V, Dina C, Hamid YH, Joly 
E, Vaillant E, Benmezroua Y, Durand E, Bakaher N, Delannoy V, Vaxillaire M, 
CookT, Dallinga-Thie GM, Jansen H, Charles MA, Clement K, Galan P, 
Hercberg S, Helbecque N, Charpentier G, Prentki M, Hansen T, Pedersen 0, 
Urrutia R, Melloul D, Froguel P: Role of transcription factor KLF11 and its 
diabetes-associated gene variants in pancreatic beta cell function. 
Proc Natl Acad Sci USA 2005, 1 02(1 3):4807-481 2. 

13. Ellenrieder V, Buck A, Harth A, Jungert K, Buchholz M, Adier G, Urrutia R, 
Gress TM: KLFl 1 mediates a critical mechanism in TGF-beta signaling 
that is inactivated by Erk-MAPK in pancreatic cancer cells. 
Gastroenterology 2004, 1 27(2):607-620. 

14. Pang YP, Kumar GA, Zhang JS, Urrutia R: Differential binding of Sin3 
interacting repressor domains to the PAH2 domain of Sin3A. 
FEB5 Lett 2003, 548(1 -3):1 08-1 1 2. 

15. Fernandez-Zapico ME, Mladek A, Ellenrieder V, Folch-Puy E, Miller L, Urrutia 
R: An mSin3A interaction domain links the transcriptional activity of 
KLFl 1 with its role in growth regulation. Embo J 2003, 22(18):4748-4758. 

16. Kaczynski JA, Conley AA, Fernandez Zapico M, Delgado SM, Zhang JS, 
Urrutia R: Functional analysis of basic transcription element (BTE)-binding 
protein (BTEB) 3 and BTEB4, a novel Spl-like protein, reveals a subfamily 
of transcriptional repressors for the BTE site of the cytochrome P4501A1 
gene promoter. Biochem J 2002, 366(Pt 3):873-882. 

1 7. Ellenrieder V, Zhang JS, Kaczynski J, Urrutia R: Signaling disrupts mSin3A 
binding to the Madl-like Sin3-interacting domain of TIEG2, an Spl-like 
repressor. Embo J 2002, 21 (10):245 1-2460. 

18. Zhang JS, Moncrieffe MC, Kaczynski J, Ellenrieder V, Prendergast FG, Urrutia 
R: A conserved alpha-helical motif mediates the interaction of Spl-like 



Calvo et al. BMC Molecular Biology 201 4, 1 5:1 0 Page 1 5 of 1 5 

http://www.bionnedcentral.conn/1471-21 99/1 5/1 0 



19. 



20. 



22. 



23. 



24. 



25. 



26. 



27. 



28. 



29. 



30. 



31. 



32. 



33. 



34. 



35. 



36. 



37. 



transcriptional repressors with the corepressor mSin3A. Mol Cell Biol 2001, 
21(15):5041-5049. 

Kaczynski J, Zhang JS, Ellenrieder V, Conley A, Duenes T, Kester H, 

van Der Burg B, Urrutia R: The Sp1-lil<e protein BTEB3 inhibits transcription 

via the basic transcription element box by interacting with mSinSA and 

HDAC-1 co-repressors and competing with Sp^.J Biol Chem 2001, 

276(39):36749-36756. 

Mathison A, Liebl A, Bharucha J, Mukhopadhyay D, Lomberk G, Shah V, 
Urrutia R: Pancreatic stellate cell models for transcriptional studies of 
desmoplasia-associated genes. Poncreotology 2010, 10(4):505-516. 
Xiong Y, Khanna S, Grzenda AL, Sarmento OF, Svingen PA, Lomberk GA, 
Urrutia RA, Faubion WA Jr: Polycomb antagonizes p300/CREB-binding 
protein-associated factor to silence F0XP3 in a Kruppel-like 
factor-dependent manner. J Biol Chem 2012, 287(41 ):34372-34385. 
Xiao TZ, Bhatia N, Urrutia R, Lomberk GA, Simpson A, Longley BJ: MAGE I 
transcription factors regulate KAP1 and KRAB domain zinc finger 
transcription factor mediated gene repression. PLoS One 201 1, 
6(8):e23747. 

Daftary GS, Zheng Y, Tabbaa ZM, Schoolmeester JK, Gada RP, Grzenda AL, 
Mathison AJ, Keeney GL, Lomberk GA, Urrutia R: A novel role of the Sp/KLF 
Transcription Factor KLF1 1 in arresting progression of endometriosis. 

PLoS One 2013, 8(3):e60165. 

Yin KJ, Fan Y, Hamblin M, Zhang J, Zhu T, Li S, Hawse JR, Subramaniam M, 
Song CZ, Urrutia R, Lin JD, Chen YE: KLF11 mediates PPARgamma 
cerebrovascular protection in ischaemic stroke. Brain 2013, 
136(Pt 4):1 274-1 287. 

Wernig M, Meissner A, Foreman R, BrambrinkT, Ku M, Hochedlinger K, 
Bernstein BE, Jaenisch R: In vitro reprogramming of fibroblasts into a 
pluripotent ES-cell-like state. Nature 2007, 448(71 5 1):3 18-324. 
CookT, Gebelein B, Belal M, Mesa K, Urrutia R: Three conserved 
transcriptional repressor domains are a defining feature of the TIEG 
subfamily of Spl-like zinc finger proteins. J Biol Chem 1999, 
274(41 ):29500-29504. 

Zhang JS, Moncrieffe MC, Kaczynski J, Ellenrieder V, Prendergast EG, Urrutia 
R: A conserved alpha-helical motif mediates the interaction of Spl-like 
transcriptional repressors with the corepressor mSin3A. Mol Cell Biol 2001, 
21(15):5041-5049. 

Bannister AJ, Zegerman P, Partridge JF, Miska EA, Thomas JO, Allshire RC, 
Kouzarides T: Selective recognition of methylated lysine 9 on histone H3 
by the HP1 chromo domain. Nature 2001, 410(6824):120-124. 
Lachner M, O'Carroll D, Rea S, Mechtler K, Jenuwein T: Methylation of 
histone H3 lysine 9 creates a binding site for HP1 proteins. Nature 2001, 
41 0(6824):1 16-120. 

Aasland R, Stewart AF: The chromo shadow domain, a second chromo 
domain in heterochromatin-binding protein 1, HP1. Nucleic Acids Res 

1995, 23(1 6):3 168-3 173. 

Chin HG, Esteve PO, Pradhan M, Benner J, Patnaik D, Carey MF, Pradhan S: 
Automethylation of G9a and its implication in wider substrate specificity 
and HP1 binding. Nucleic Acids Res 2007, 35(21):731 3-7323. 
Fuks F, Hurd PJ, Deplus R, Kouzarides T: The DNA methyltransferases 
associate with HP1 and the SUV39H1 histone methyltransferase. 
Nucleic Acids Res 2003, 31 (9):2305-231 2. 

McConnell BB, Yang VW: Mammalian Kruppel-like factors in health and 
diseases. Physiol Rev 201 0, 90(4):1 337-1 381 . 
Havula E, Teesalu M, Hyotylainen T, Seppala H, Hasygar K, Auvinen P, 
Oresic M, Sandmann T, Hietakangas V: Mondo/ChREBP-Mlx-regulated 
transcriptional network is essential for dietary sugar tolerance in 
Drosophila. PLoS Genet 2013, 9(4):e 1003438. 
Small KS, Hedman AK, Grundberg E, Nica AC, Thorleifsson G, Kong A, 
Thorsteindottir U, Shin SY, Richards HB, Consortium G, Soranzo N, Ahmadi 
KR, Lindgren CM, Stefansson K, Dermitzakis ET, Deloukas P, SpectorTD, 
McCarthy Ml, Mu TC: Identification of an imprinted master trans regulator 
at the KLF14 locus related to multiple metabolic phenotypes. Nat Genet 
2011,43(6):561-564. 

Gebelein B, Fernandez-Zapico M, Imoto M, Urrutia R: KRAB-independent 
suppression of neoplastic cell growth by the novel zinc finger 
transcription factor KS^. J Clin Invest 1998, 102(1 1):1 91 1-1919. 
Gebelein B, Urrutia R: Sequence-specific transcriptional repression by KS1, 
a multiple-zinc-finger-Kruppel-associated box protein. Mol Cell Biol 2001, 
21(3):928-939. 



38. Seo S, Lomberk G, Mathison A, Buttar N, Podratz J, Calvo E, lovanna J, 
Brimijoin S, Windebank A, Urrutia R: Kruppel-like factor 11 differentially 
couples to Histone Acetyltransferase and Histone Methyltransferase 
Chromatin remodeling pathways to transcriptionally regulate Dopamine 
D2 receptor in neuronal cells. J Biol Chem 2012, 287(1 6):1 2723-1 2735. 

39. Fernandez-Zapico ME, Lomberk GA, Tsuji S, DeMars CJ, Bardsley MR, Lin YH, 
Almada LL, Han JJ, Mukhopadhyay D, Ordog T, et al: A functional 
family-wide screening of SP/KLF proteins identifies a subset of suppressors 
of KRAS-mediated cell growth. Biochem J 201 1, 435(2):529-537. 

40. Lomberk G, Bensi D, Fernandez-Zapico ME, Urrutia R: Evidence for the 
existence of an HP1 -mediated subcode within the histone code. 
Nat Cell Biol 2006, 8(4):407-415. 

41. Lomberk G, Mathison AJ, Grzenda A, Seo S, DeMars CJ, Rizvi S, 
Bonilla-Velez J, Calvo E, Fernandez-Zapico ME, lovanna J, Buttar NS, Urrutia 
R: Sequence-specific recruitment of Heterochromatin protein 1 via 
interaction with kruppel-like factor 11, a human transcription factor 
involved in tumor suppression and metabolic diseases. J Biol Chem 2012, 
287(1 6):1 3026-1 3039. 

42. Irizarry RA, Hobbs B, Collin F, Beazer-Barclay YD, Antonellis KJ, Scherf U, 
Speed TP: Exploration, normalization, and summaries of high density 
oligonucleotide array probe level data. Biostatistics 2003, 4(2):249-264. 

43. Prasanth SG, Prasanth KV, Siddiqui K, Spector DL, Stillman B: Human Orc2 
localizes to centrosomes, centromeres and heterochromatin during 
chromosome inheritance. Embo J 2004, 23(1 3):265 1-2663. 

44. Kondo Y, Shen L, Ahmed S, Boumber Y, Sekido Y, Haddad BR, Issa JP: 
Downregulation of histone H3 lysine 9 methyltransferase G9a induces 
centrosome disruption and chromosome instability in cancer cells. 
PLoS One 2008, 3(4):e2037. 

45. Wettenhall JM, Simpson KM, Satterley K, Smyth GK: affylmGUI: a graphical 
user interface for linear modeling of single channel microarray data. 
Bioinformatics 2006, 22(7):897-899. 



doiilO.1 186/1471-2199-15-10 

Cite this article as: Calvo et al.: Single and combinatorial chromatin 
coupling events underlies the function of transcript factor kruppel-like 
factor 1 1 in the regulation of gene networks. BMC Molecular Biology 
2014 15:10. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



o 



BioMed Central 



